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Abstract: In this paper, we prove the existence and uniqueness of local strong 
solutions of the hydrodynamics of nematic liquid crystals system under the initial 
data satisfying a natural compatibility condition. Also the global strong solutions 
of the system with small initial data are obtained. 
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1 Introduction 

It is well known liquid crystals are states of matter which are capable of flow and in 
which the molecular arrangement gives rise to a preferred direction. By Ericksen-Leslie 
theory, the compressible nematic liquid crystals reads the following system: 



pt + div (pu) = 0, 

1 r, 

{pu)t + div(pn (^u) + s/p = pAu - Adiv(vrf V^^ - o (I V^^l + F{d))I) 
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in r2 X (0,T), for a bounded smooth domain Q C M'^. p{t,x) is density, u{t,x) the 
velocity field, x) orientation parameter of the liquid crystal and p{p) pressure with 

p = p{-)eC'[0,oo), p(0) = 0. (1.4) 

The viscosity coefficients /i. A, u are positive constants. The unusual term \/d \/d 
denotes the 3x3 matrix whose {i,j)-th element is given by Y^k=i^^i^k^^i'^k- ^ 
unite matrix. 

f{d) = l(Mr - l)d and F{d) = ^(l^r - 1)'. 
We are interested in the initial data 

p(0, x) = po > 0, u(0,x) = wo, (i(0, x) = (io(a;), Vx G (1-5) 
and the boundary condition 

w(t,x) = 0, d{t,x) = do{x), |cio(2;)| = l, V(t, x) G (0, T) x (1.6) 

The first author pQ firstly considers the model fll.ip - fll.4l) and gives a global existence 
of finite energy weak solutions by using Lions's technique (see [2||3] ). The incompressible 
model of a similar simplified Ericksen-Leslie model has been studied in the papers [IHH]. 
Recently, for the incompressible Ericksen-Leslie model, Wen-Ding [5] have gave the 
proof of local strong solutions in two and three dimensions, and Lin-Lin- Wang [10] 
have established the existence of global (in time) weak solutions on a bounded smooth 
domain in two dimensions . 

Throughout this paper, we use the following simplified notations 

In order to obtain strong solutions, we need {po,Uo,dQ) satisfies the regularity 

po G Mo G Hl^ n H\ do G H\ (1.8) 
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and a compatibility condition 

fiA uo- Adiv(vrfo Vdo - ^{\vdof + F{do))I) - Vpo = Vx G Q, (1.9) 
for some g E LF' . 

The compatibility condition is a natural request to insure ||y/p?/t(0)||L2 bounded, and 
a compensation to the lack of a positive lower bound of po- It is firstly introduced by 
Salvi, Straskraba |llj and the authors of [12] independently. 
Our main result is the following theorem: 

Theorem 1. Part I: (Local existence) Under the assumptions of the regularity 
condition (11. 8p and the compatibility condition (II. 9p . there exists a small T* G (0,T) 
such that the system (ll.ip - (ll.6p has a unique local strong solution {p,u,d) satisfying 

peC{[0,T*]-W^'^), pteC{[0,T*];L^), 

ueC{[0,T*];H^r]H^)r]L\0,T*;W^^^), d e C {[0 , T*]; H'^) , 

ut G L2(0, T*; Hi), dt G C([0, T*]; n L\<d, T; H^), 

^ut G C([0,T*];L2), du G L\0,T*-L''). 

Part II: (Continuity of initial data) Suppose that {p,u,d) is another solution of 
(ll.ip -( n~4p with the following initial boundary conditions 

p{0,x) = pq>0, u{0,x) = uo, d{0 , x) = doi^x) , WxeQ, 
u{t, x) = 0, d{t, x) = do{x), V(t, x) G (0, T) x dQ, 

then for any t G (0,T*], the following quantities tend to zero when (j)Q,Uo,do) — )■ 
{po,Uo,do) in W^'^ x x H^: 

\\d - d\\H^{t), \\d - d\\L2{o,T*;H3), \\ dt - dt\\ (t) , \\dt - dt\\L^(o,T*;H'>^)' 
\\U-U\\m{t), \\u-u\\l2(0,T*;H^), \\y/p{Ut- Ut)\\L^{0,T';L2), || P - p|| L6 (t) . 

In general, as Navier-Stokes equations we have not the global existence with large 
initial data. However, a blow-up criterion is obtained in our article [13] and for small 
initial data, we have the following global existence. 
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Theorem 2. (Global existence with small initial data) Let a be a nonnegative 
constant and m a constant unite vector in M^. Then there exists a positive constant 
small such that if the initial data satisfies 

max{||po - o;\\w^,6, ||no||j/2, ||cio - ^11^3, WgWh} < 0, (1.11) 

for all 9 G (0, 9q], then the system fll.ip - fll.6l) has a unique global strong solution. 

Remark 1. In this paper, we only consider the case a = 0, because the other case 
{a > 0) can be induced to the problem with a positive initial data of density. 

The methods to prove the Theorem [T] and Theorem [2] are the successive approxima- 
tion in Sobolev spaces (see [1H[I5]). We generalize this method to two variables, which 
is based on the following careful observation on coupling terms of d and u: 

(1) ||(i||^2 can be deduced by ||n||^2dr and 

(2) ||'u||^2 can be derived by HdH^adr. 

Since the initial density has vacuum, the movement equation (11.21) becomes a degenerate 
parabolic-elliptic couples system. To overcome this difficulty, as usual, the technique is 
to approximate the nonnegative initial data of density by a positive initial data. For a 
special linear equations of the system (ll.ip - (ll.3l) . we prove not only the local existence 
of strong solutions with large initial data but also the global existence of strong solu- 
tions with small initial data. Employing energy law and higher energy inequalities, we 
can prove both the uniqueness and the continuous dependence on the initial data. 

This paper is written as follows. In section [2], after establishing a linear problem 
of the nonlinear problem (ll.ip - (ll.3l) . we prove local existence of a strong solution to 
the special linear problem with a positive initial data of density. And we also establish 
some uniform a prior estimates, which imply the existence of a local strong solution to 
the linear problem when the initial data of density allows vacuum in section 12.31 In 
section [31 after constructing a sequence of approximate solutions, a strong solution of the 
nonlinear problem (ll.ll) - (ll.3p is obtained. In section HJ the uniqueness and continuous 
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dependence on initial data are proved. 

2 A linear problem 

2.1 Linearization 

At the beginning, we linearize the equations fll.ip - fll.3p as following: 

pt + div(p^;) = 0, (2.1) 

dt + v- \/d = u{Ad ^{n + m) ■ {d — m)n), (2.2) 

put + pvu + \7p{p) = pAu- X{\ydf{Ad - f{d)), (2.3) 

where m is a constant unite vector and {v,n) satisfies the following regularity 

V e C([o, T]; n n l\o, t- W^^), vt e L\o, T- Hi), 

n e C([0, T]- H^) n L2(0, T- H^), nt G C([0, T]; H^) n L\0, T; H^). 
Assume further that 



(2.4) 



(2.5) 



t'(0, x) = Mo(a;), n(0, x) = (io(2^)) Vx G fi, 

v{t, x) = 0, n{t, x) = do{x), V(t, x) G (0, T) x dn. 

2.2 The existence of approximate solutions 

For each 5 G (0, 1), let be the solution to the boundary value problem: 

pA4- Adiv(v(/o Vdo - Vt^ol' + F{do))I) - p{pi) = {pf,)h, 
ul = 0,Vx G dQ, 

where Po = Po + Then — uq in fl as 5 — t- 0. 

For the linear equations f l2.ip - fl2.3p . we have the following theorem. 

Theorem 2.1. Assume po, Uq and do satisfy the regularity U.8\) and {v,n) satisfies the 
above conditions ^2.4^ and Ii2.5\) . Then there exists a unique strong solution {p,u,d) 



(2.6) 
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of the linear equations Ii2. l\) - ^27^) with the initial data (pg,Mg,(io) one? the boundary 
condition U.6\) such that 

peCi[0,T];W'''), p,eC([0,T];L6), 

u e C([0, T]; Hi n H') n L2(0, T; t^2,6)^ „^ ^ c'([0, T]; L') n L2(0, T; //g^), 
G C([0, T]-H^) n L2(0, T; Z/^)^ G C([0, T]; i/i) n ^^(O, T; if^)^ 

gL2(0,T;/7-1), dueL\<d,T-L^). 

We will prove the theorem by the following three lemmas. Firstly, suppose the 
constants cq, ci and C2 satisfying 

Co > 1 + IIpoIUm + Nolk^ + Ihob^ + lkl|L2, (2.8) 
ci > sup + ||n||^i + ||ni||j;^i) + / (|| V lli^ + ll^ll^^.e 

0<t<T " Jo 

+ \\V^rit\\l^ + M\H-Mt, (2.9) 
02 > sup (II v^^I|l2 + II V^'^IUO, (2-10) 

0<i<T 

C2 > Ci > Co > 1. (2-11) 

Lemma 2.1. Assume po and v satisfy the regularities ( [j.<§)) anc? ( [^.^[ j respectively. Then 
the problem Ii2.1\) and U.5\) has a global unique strong solution such that 

peC{[0,T]-W''^), aGC([0,T];L6). (2.12) 
Moreover, if v,n satisfy I{2.8\) - I[2J1\) . then there is a Ti = min{c^\T} such that \/t G 



kM(t)<Cco, ||pt|U6(t)<CcoC2. (2.13) 
In particular, 

IbllLs(t) < CAf(co), II VJ5|U6(t)<CM(co)co, |bt||L6(t) < CAf(co)coC2, (2.14) 
where the constant M{cq) is defined by 1^2. 18\) . 



Remark 2. The above lemma is also true to the problem (12.11) with the initial data 
and the estimates 12.13112.1^ also hold for all small 5. 

Proof. Now let's start to prove Lemma 12. 1[ The existence is obvious, based on the 
classical method of characteristics (see [15]). So 

p{t,x) = po{y{0,t,x))exp{- divt;(r, ?/(r, t, x))dr), (2.15) 

^0 

where y{T,t,x){E C{[0,T] x [0,T] x Q)) is a solution of the initial value problem: 

-§;:y{T,t,x) = v{T,y{T,t,x)), < r < T, 
y{t, t,x) = X, < t < T, X E Q. 

Moreover, we have 

\\p\\w^,6{t) < \\po\\wifiexp{ \\\/ v\\wi.e{T)dT), VtG[0,T]. 
So that for all t e [0,Ti], 

IIpIIm/i.sW < ||Po|Ui.6 exp( / W \/ v\\wh6{T)dT) 







Because of ||p||l°° < C||p||w^i'6; set 



We can obtain for all t G [0, Ti], 



\pt\\Loit) = \\ - V ■ \/p - pdivvWie 

< \\v\\l^\\ V pIIlc + V 

< Il^^ll/f2||p|ki.6 + llplki.Hl^lk^ 

< CcoC2. 



(2.16) 



<Coexp(t2(/ II V ^^llH/i>6(r)dr)2) 

^0 

< eco. (2.17) 



M(co)= sup (l+p(-)+p'(-))- (2.18) 

0<-<Ceco 
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Similarly, for all t G [0,Ti], 

MAt) < CM{co), II VPWAt) < CM(co)co, 

WPtWAt) < CM(co)c2Co. 

□ 

Lemma 2.2. Under the hypotheses of Theorem \2.1[ then the equation ^2.2^ with the 
initial boundary conditions / fi.5]) -/ l77^) has a global unique strong solution d satisfying 



Moreover, ifv,n satisfy I[2.8\) - I[K11\) . then there is a = minjcg ,T} such that 
sup (||d||Hi +crl V'^ dh^ + 4\\dt\\m + c^'^c^^W V d\\H^) 

0<t<T3 

|2 , Jii j||2 1. ^ nJ- 



+ / \\dtrm + cm\H-^^t<Cci (2.19) 

Proof. By the classic Galerkin method to the linear parabolic equation (12. 2 p with (II. 5p - 
(II. 6p . the existence and regularity of d described in (12. 7p can be obtained. 

Differentiating fl2.2l) with respect to time, multiplying by dt and then integrating 
over f2, we can deduce that 

--T- /" MiPdx + z/ [ Iv^tPdx 

<C'(lktlU«ll y d\\L2\\dt\\L-i + ||^^||l°°|| y dt\\LA\dt\\L^ + ||^||L3||nJi6||(i-m||i3||djL6 
+ ||n + m||L3||nt||j:6||(i - m||j:3 ||(it ||l6 + ||n||j:6||n + m||i6||dt||^3) 

5 

= Y.h. (2.20) 



i=l 

Here 



h < Cri\\vt\\U\ V d\\l. + Cr/"i ll^/^lli, + ^|| y «?t||i., 
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<C||t;||?,.|M,||i. + ^|| V^dli., 



/3 + /4 < C((||n||i. + l)(||n||^. + l)||nJia||d - m||i. + || y ^lli. + \\d- m||i.) 



+^\\ydt\\i2, 



h < C||n||^i||n + m||^i||dt|||2 + ^|| y d,||2,^ 



where the small positive constant rj will be fixed later. 
Since 

/ \ \/d\^dx = 2 I \/d:\/dtdx< [ \\/d\'^dx+ [ \\/dt\^dx, (2.21) 

[ \d-m\'^dx< [ \d-m\'^dx+ [ \dt\'^dx, (2.22) 
Jn Jn Jn 

combining flZ^ . (ES]) and (ESSD, we get 

^ ^ ' ' + \\/ dl"^ + \d-m\'^dx+ I \\/ dt\'^dx 



/ Mtr 



|2 

1^1 



dt 

<c(ll V ^^11^2 + - Hh + IMtlli^) ■ (^"^ + + ll^ll?/! 11^ + ^\ 

+ + (ll^lli^ + mn\\]j. + 1)11 V ^illt^ + !)• (2.23) 

From the equation (12.21) . we can deduce 



kWh^^) <C(|| A dolU^ + IImoIIhHI V doWh^ + No + "^lUHIt^o - "^IIl^ llt^olUO 

<Cc|j. (2.24) 

Hence, by Gronwall's inequality, we can deduce from (I2.23P 

\dt\'^ + \\/ d\^ + \d-m\^dx+ j j \\j dt^^dxdr 



< Cclexp{C [ V ^ + \Mh2 + Wnfjji + \\n + mfjj^ + ri\\vt\\Hi 
Jo 

+ i\\n\\h + mnrH. + mVnt\\hdT). (2.25) 



Taking rj = c^^ and using the assumption (I2.8p - (l2.1ip . we obtain 



sup [ \dt\'^ + d\'^ + \d-m\'^dx+ [ [ \ dtl'^dxdr < Cc^, (2.26) 
o<t<T2 Jn Jo Jn 

where T2 = min{c^®, Ti}. 

From the equation (12. 2p and using the elliptic estimates, we get 

\\d - m\\H-2 < C{\\do - m\\H2 + WdtWi^ + \\v ■ V^^IIl^ + \\{n + m) ■ {d - m)n\\L2) 

1 1 

< C(cq + II V II //I II V '^IIl2||c? - ^Wh^ + \\^ + "^IIlsIM ~ "^IIlsII'^IUO 

< C(co + ||v||hi|| V c^IIl2 + ||n + m||j^i||d- m||j^i||n||j^i) + ^|M||/^2. 
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So 

II V' d\\L2 < C{cl + clcl + clcl) < Cclcl (2.27) 

Applying the operator y to the hnear equation fl2.2p . we get 

p A {\/d) = \/dt + v(^' ■ Vc?) + 4 V (('^ + m) ■ (d - m)n). (2.28) 

By the elhptic estimates, we can estimate the term || v ('^ ~ ^) 11/^2 as follows 

II V {d-m)\\H2 <C(|| y dtWi^ + II V ■ Vd)\\L^ + 11^ V {{n + m) ■ {d-m)n]\\L2 
+ No - mWu-i) 

- - 

<C{\\ y dt\\L2 + II yv\\HA\ y d\\l,\\ V chilli + hWuAl y d\\L2 

+ II V ?T-||//i||?^lki||rf- "^||//i + \\n\\]i2\\ V d\\L2 + \\dQ - m||//3) 
<C{\\ydA\L- + clcicl), (2.29) 

where we use the assumption fl2:8|) - (l2ni|) and (12:261) - dXTTD . 

Differentiating (12. 2p with respect to time and taking inner product with Arf^, then 
we can derive 



Id/", ^ ,2, 



< ^— / \ \j dt\^dx + u I \Adt\^dx 



2dt 



/ (tJf ■ v^) ^ f^tda: + {v ■ \jdt) A dfdx H — ^ / {rit ■ {d — m))n A dtdx 
Jn Jq ^ Jn 

+ — / (n ■ dt)n A (i^dx H — - {{n + m) ■ {d — m))nt A dtdx 



5 

= E"^^' (2.30) 
where we use the elliptic estimate ||(it||^2 < C'H A dtjl^a- 
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Here 

|Ji| = \ {vf \7d) A dtdx\ 
Jn 

< II V VtWi^W V o^IIl«|| V dth^ + \\vt\\LA\ d\\L6\\ V dt\\L2 

< ^11 V ^tlli^ + V^'^W V dWljiW V dth^W V f^tll^i + ^Iktllisll V dt\\l2 

< V\\ V ^tlli^ + Cv-^e-'clcm V ^^illi^ + 4dt\\h + v\\vt\\m\\ V dt\\l2 

+cv-\\\vdt\\h + 4cl4), 

IJ2I = \ {v-\7dt)Adtdx\<Ce-'^cl\\\ydt\\l2 + e\\Adt\\l2, 
Jn 

\M < Ce-in||i^||(i,||i.+£||ArfJi. <Ce-ic^c^ + e|| AdJi., 
I^sl + I^sl < C6-\\\n\\l^ + l)\\nt\\l2\\d-m\\l^+e\\Adt\\l2 
< Cs-^cl4cl + s\\ Adt\\l2, 

where the small positive constants e and rj will be fixed later. 

Substituting the above Ji — J5 into fl2.30p and taking e small enough, we get 

- / \x7dt\Mx+\\dt\\l2<Ar,{t)\\x7dt\\l2 + B,{t), (2.31) 

where 

A, it) = Civ-'4cl' + v\\vt\\l.+V-' + cl), 

B, it) = C(7^||t;,||^. +ri4c?c^ + r'c^c?c^ + 4c^). 

In view of the inequalities fl2.26p - p.27p and taking rj = Cg^, we have for all t G [0,T2], 
/o Ar,{s)ds <C + CclcfclH, /o B^{s)ds <C + Cclc\clt. (2.32) 
From the equation fl2.28p . we can estimate the term || v '^t('^)||L2 at time 

II V c^t(o)||L2 < C'dl^o - "^1^3 + II V ^oIIlsII o^o||l6 + II^^oIIl^II ^^olUe 

+ 11 V 4||L6||f^0 - "^lUsll^olUe + IMolliejl V C^olUe) 

< Cel. (2.33) 
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Applying Gronwall's inequality to f l2.3ip and using fl2.32p -( l2l33ll . we can obtain 

.T3 



sup /" I V dtl^dx + [ ' WdtWl^dt <C{cl + 444^) exp {C + Ccjclcl^n) 
o<t<T3 Jn Jo 

<Ccl (2.34) 

From (I2:29|) and ( 1234]) . we have 

II V d\\m{t) < C(c| + clci4), \/t e [0,n]. (2.35) 

Hence 

ft 

j,3dT<C, \/te[0,Ts]. (2.36) 



Since ([OS]), dOZD, (EMD and (lO^ can deduce the estimate (1^19]) . we complete the 
proof. □ 

Lemma 2.3. Under the hypotheses of Theorem \2.1\ then there exists a global unique 
solution u of the equation /{2.3\) with the initial data Uq and the boundary condition 
( fi.61) . and the solution u has the regularity in {2.1). 



Moreover, if v,n satisfy li2.8\) - f2Jl\) . then, for all small 6, 

sup (M(co)4'||m||j^i + M{co)4cf\\ v' u\\l2 + M(co)4i VP^tlU^) 

0<t<T3 

+ r ^oll V UtWh + \\u\\l.,,edt < C4'M\co). (2.37) 
Jo 

Proof. Because Po > S > 0, it follows from the representation fl2.15p that 

p{t,x)>Sexp{- [ \\7 v{T)\w^,edT) > S, V(t, x) G [0, T] X (2.38) 
Jo 

where 5 is a positive constant. 

Thanks to fl2.38p . we change (12.30 into the following form: 

ut + v-^u+-^p=-Au - -{\ydf{Ad - f{d)). 
P P P 
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Applying the Galerkin method again to the above equation with the initial data Uq and 



the boundary condition ( ll.6p . we can deduce the existence and regularity of u described 
in (ETD. 

Differentiating (12. 3p with respect to time t, multiplying by Ut and then integrating 
over Q, we can derive 

= {- V Pt - Ptv ■ - ^pv ■ \/Ut - pvt ■ \/u))utdx 
Jn 

-A / {\7dtf{Ad-f{d))utdx-X [ {xjdf{/\d- f{d))tutdx 
Jn Jn 

6 

= (2.39) 

k=l 

Here 

\Ki\ < C\\p,\\l. + ||| V < CM\co)clcl + ||| V ^tlli., 

\K2\ < C\\pt\\l4v\\l4 V ^^11^2 + ||| V utWl^ < CclclclW V u\\l, + |ll V Ml^. 

\K,\ < C\\p\\L^\\v\\U\^pu,\\l, + ^11 V u,\\l. < CclcW^pmWl. + ^11 V n,\\l,, 
\Ki\ < v\\ y Vt\\\2\\^ut\\l2 + r]'^\\p\\LA\ V^^IIlHI yAm 



< v\\ V vt\\i2\\y/put\\i2 + V c^ll V u\\i2 + 



li^sl < C(|| v' 4^4 V dt\\L4ut\y + II V d4L4d\\L4ut\\L4d\\lo. 

+ 11 V (itllLalMIUHI^tlUs) 

< c\\ d\\,4 rflUall V dtWh + C\\ V d,\\l2\\dU4d\\H4drH2 
+C\\vdt\\l2\\dh4d\\m + ^\\VUt\\l. 

< Cclclcf + f^\\s7u,\\l2, 
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l^el <C(|| d\\L4 V dt\\L4ut\\LB + II V ^^lU-ll V ^^tllL^H V uth^ 

+ II V d\\L^\\dt\\L3\\Ut\\L^ + II V 0?|U3||dt||i6||rf||^6||tXt||i6) 

<c(|| c^IIlHI ^^IIlbII V dtWh + II V t^ll^i,6|| V f^tlli^ 

+ \\dt\\HA\dt\\L4v d\\l2 + II vo^IIlHI Vc?IIhi|| Vc^illi^ 



II ii2 

+ ^llv^.lL. 



On the other hand 
d 



^11 V* 



/ S7U : S7Utdx < C\\ V u\\l, + ^11 V UtWl^. (2-40) 
Jn ^ 

Substituting li^il — \Kq\ into fl2.39p and combing with (12.401) . we obtain 

[ {p\ut\'^ + \ \y u\'^)dx + [ \\7ut\'^dx 
'^t Jn Jn 

< ^.mW^utWh + II V u\\l) + ^(t) + C|| V ^ll?.i, (2.41) 

where 

^r,{t) = C{4clcl + v\\ V ^tlli. + V-'cl), .m = C4cicfM\co). 
Taking r] = c^^ and using the estimates fl2A3D - ( 1271^1) and fl219l) . we have Wt e [0, Tg], 

£^^{s)ds <C + Celt, /J ^{s)ds < Cc\c\cfM^{co)t + C. (2.42) 

Multiplying f l2.3p by Mj, then integrating it over Vt and using the Young's inequality, we 
can obtain 

/ p\ut\'^dx{T) < C p\v\'^\\/u\'^ +\p~^fiAu- Xdw{\/dQ\/d 
Jn Jn 

-li\vd\' + Fid))I)-Vp)\'dxiT). 



Hence 



hmsup / p\u,\MxiT) < Cicl + \\g\\l,) < Cel. (2.43) 

T^0+ J 
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Integrating f l2.4ip with respect to time (t, t) and letting r — )■ 0+, thanks to fl2.43p and 
Gronwall's inequahty and the inequahty fl2.42p . we have for all t G [OjTa], 



{p\ut\' + 



JoJq 



< CclM\co) + C [ II v^ill^idr. 
Jo 



(2.44) 



Using the elliptic regularity result to the linear movement equation f l2.3p . we can 
estimate the term || "^lU^ as follows 

Wv^Wrn < C{\\put\\L2 + \\pv\yu\\L2 + \\\7p\\L2 + \\\7u\\L^ 

+ ||Mr(Arf-/(d))|UO. (2.45) 

From the assumption (EH]) -([2111]) and the estimates fl213l) . fl2:i4l) and fl2J9l) . we can 
derive 

Wi^df A d\\L2 < C\\{s^df{dt + V ■ yd)\\L2 + C\\{ydf[{n + m)-{d- m)]n\\L2 

< V dWl^W V d\\jji\\ V dth^ + II V o^ll//i||^^|ki|| V d\\m 
+ ||(i||j:^2||77, + m\\L6\\d- m||//2 Unlike) 

< 044, (2-46) 

UvdffmL^ < C\\x7d\\H4\\d + m\\L4d-m\\L4d\\H^<C44^ (2.47) 
and 

Wputh'^ + IIp^^ • v^lU^ + II vpIIl^ + II V ^IIl2 

< \\Vp\\L^\\^/pUt\\L2 + ||p||L-||f IIlsII V + CM{cq)Co + \\ S/ u\\l2 

1 11 

^IIl^II V ''^IIhi + cm{cq)co + II V "^lU^ 

< CM{co)co + C4{\\s7u\\L^ + \\^ut\\L^) + ^\\Vu\\m- (2-48) 
Taking (|226])-(|228]) into ([235]), we can deduce 

II V < €44' + CM{co)co + C4{\\ V u\\l'2 + W^ikh^)- (2.49) 
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Taking (12 .49^ into (12 .44^ and using Gronwall's inequality, we can derive for all t € [0, T3] 

[ {p\ut\^ + \\yu\^)dx+ [ [ \\yut\MxdT<Ccl°M\co). (2.50) 
Jq JoJq 

So substituting (1230|) into fICTD . we get 

II V^lli/i <Cc?4°M(co). (2.51) 



Using the elliptic regularity result to the linear movement equation (l2.3p . the term || 
^11^6 can be controlled as follows 

II V^^IIls < C{\\put\\L6 + \\pv ■ v^IIls + II \7p\\l6 + \\{\7df{Ad-f{d))\\Le 

+ ||Vw||lO- (2.52) 

It follows from the assumption fra- fimi) and the estimates fl^TT^ . fl^THD and fICT]) 
that 

llpMtlUe + \\PV ■ V^IIls + II V Pha + II V uWlo 

< ||p||L->||Mt||L6 + ||p||oo||v||l°-|| V + CM{co)co + CCiCo°M(co) 

< Cco\\s^Ut\\L^ + Cc2clcl^M{co). (2.53) 

It follows from the estimates (12.271) and (I2.35P that 

\\{S7df{Ad-f{dmLe < C||rf||^3+C|| v^^l|L6||d-m||H2|M + m||H2|M||^^2 

< C4c^^cl\ (2.54) 

Taking (I2.53P and (12.541) into (I2.52p . integrating it over time and using the estimate 
(12.501) . we can derive for all t E [0, T3], 

/ II v' "lliadr < Ccl [ II V utWl^dr + C^^M\co) < C^^M\co). 
Jo Jo 

□ 

It is obvious that Lemma I2.1j -Lemma 12.31 imply Theorem 12.11 
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2.3 Local existence of a solution to the linear problem (12.11 )- 



(123]) 



Since the estimates obtained in the Lemma l2.H -Lemma 12.31 are uniform for all small 
6, we have the following theorem: 

Theorem 2.2. If the initial condition {po,Uo,do) satisfies the regularity i\1.8^ and the 
compatibility condition ( fi.gj) . then there exists a unique strong solution (p, u, d) to the 
linear equations with initial boundary value U.5\) - I[l~^) such that 

ueC{[0,Ts];H',nH^)nL\0,Ts;W^'^), Ut e L\0,Ts; H^), ^ ^ 

(2.55) 

d G C([0, T3]; H') n L2(0, Tg; H^), dt E C([0, T,]; H',) n ^^(0, T,; H'), 

Moreover, {p,u,d) also satisfies the inequalities /12.1^)-(S.14\)> ^KIW and [KY^. 



Before proof, we give the following two classical lemmas which are proved in the 
book m. 



Lemma 2.4. Let Y = {v\v G L"«(0, T; Xq), G L"i(0,r;Xi)} with norm \v\y = 
k|L'^o(o,T;jfo) + kt|L°i(o,T;Xi) whcrc 1 < Oq, «! < oo,Xo C X C Xi arc Banach spaccs and 
Xo,Xi are reflexive. Suppose that the injections Xq X Xi are continuous, and 
the injection from Xq into X is compact. Then the injection from Y into L"°(0, T; X) 
is compact. 

Lemma 2.5. (An Interpolation Theorem) Let V, H, V be three Hilbert spaces, each 
space included in the following one as V G H = H' G V , V being the dual ofV. If a 
function u belongs to L^(0, T; V) and its derivative Ut belongs to L'^{0, T; V), then u is 
almost everywhere equal to a function continuous form [0,T] into H. 

Proof. We now start to prove Theorem 2.2. From the Theorem 12.11 and the estimates 
(12.131) . (I2.14p . (I2.19P and (I2.37p . by the compactness Lemma YTM. there exists {p,d,u) 
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such that 

d\ u^) {p, d, u) in L2(0, Tg; U x x W^'''), (2.56) 
/ 4 p in L°°{0,T3;W^'^), as 5 -> 0, (2.57) 

where Vr G (1, +oo) and Wq G [2, +oo). Hence p = p, a.e.. 

Because of the lower semi-continuity of various norms, the estimates f l2.13p . f l2.19p 
and (12.37P also hold for {p,u,d). So for almost every (t,x)(G [OjT^] x Q), {p,u,d) 
satisfies the system (12. ip - ( 1213]) which means (p, u, d) is a strong solution to the linear 
equation ( I2.ip - (]2.3p with initial-boundary conditions ( ll.5p -( !L6l) . 

The solution (p, u,d) is unique: From the Lemma 12. H p is the unique solution of 
the linear equation (12. ip . Using the same method as section HI we can prove d and u 
are the unique solution to the linear equations (12. 2 p and ( 12. 3p respectively. 

Finally, we will prove the time continuity of the solution (p, d). The solution from 
the Lemma 12.11 is the same as from the approximation (I2.56P due to the uniqueness of 
solution. So we get 

pG C([0,T];H^^'6). (2.58) 

From the linear equation (12. ip . we easily show 

PiGC([0,T3];L6). (2.59) 

By the interpolation Lemma 12.51 (I2.37P can deduce that 

dt G L\0,n;H^), d G L^{0,Ts;H') ^ d e C{[0,n]; H^). (2.60) 

Differentiating the linear equation ( 12. 2p with respect to time and space, we get 

Vdtt + V(^ ■ Vd)t = viS7 ^dt- j2\/[{.n + m) -{d- m)n]t). 

By the estimate (I2J9|) . we deduce ydu G ^^(0, T3; R-^). Because \/dt G ^^(0, Tg; H^), 
by the interpolation Lemma 12.51 again, we have 

rf< GC([0,T3];i/o'). (2.61) 
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From the linear equation fl2.2p . we get 

AdeCi[0,n];H'). (2.62) 
By the interpolation Lemma [2.5[ f l2.37p can deduce that 

ut e L\0,Ts;H'), u G ^^(O, Tg; W^^.e^ ^ ^ ^ C([0, Tg]; //q')- (2-63) 
From the linear equation (12 .Sp and the estimates (I2.13p . (12.190 and (I2.37p . we obtain 



{putApUt)t) e L\0,T,-H') X L\0,T,-H-') put e C([0, Tg]; L^). (2.64) 

From the linear equation (12. 3p and the elliptic regularity estimate | \/'^u\l2 < C\ Au\l2, 
we have 

ueC{[0,n];H^). (2.65) 
So we get time-continuity of {p,u, d) from (I2.58p - (l2.65p . □ 



3 Iteration and existence in Theorem 1 

Set 

c, = Ccl'M\c,), C2 = cl T, = min{c^^\T}. (3.1) 

At the beginning, let's choose a initial data of iteration {vP{t,x),d^{t,x)). vP{t,x) 
satisfies the following heat equation 

0t - A0 = 0, with (j)\t=Q = uq, (j)\on = 0, 

and d^it^x) = do{x). Because {ci,C2,Ts) depends only on cq, we can choose a T*(g 
[0, T3]) so small that fl^ - fim]) hold for m° and d° with instead of T. 

Replacing {v,n) by {u^,d^) and using the Theorem 12. 2[ we can obtain the solution 
{p\u\d^) of ([2IID-([23D with ([L5])-([L6D, and it satisfies the estimates (l27[3D . (ICTj) 
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and f l2.37p . Inductively, for all k G A^^, replacing {v,n) by {u^~^,d''^^) and using the 
Theorem 12.21 we can obtain a sequence {p^,u'',d^) of solution of fl2.ip -f l2l3|) with (11 .Sp - 
i ll. 6p . and they satisfy the following estimates with the same (cq, ci, C2, T^,) independent 
of A; e N+: 

sup (ii^^ii^. + \\d'y^ + + c^%\\ u'w,. + II rflUO) 

o<t<n 

+ [ II V ^'lli^ + \\u'\\w2,, + II V' d^Wh + Wd'Wlsdt 
Jo 

< CcfM\co) (3.2) 

and 

sup (||p^|Ui,6 + IIp^IUs) < Cc2Co, sup llv^M^ ||l2 < CM{co)cl, (3.3) 
o<t<r* o<t<T, 

sup (|b1ki>6 + Ibt IIlO < CM(co)c2Co, sup II V ^^1^2 < Cclcicl (3.4) 
o<t<T* o<t<r. 

We will show [p^^u^^d^) converges to a strong solution to the original nonlinear 
problem ffLT|) -fOD. 
Define 

-k+i ^ pk+i _ pk ^ f+^^^k+i_^k^ -k+i ^^k+i _^k_ (3.5) 

Since {p'^.u^.d'^) and {p^^^ .v!^^^ ,d'^^^) satisfy the linear equations (12. ip - (12. 30 . we 
have 

Pi^+^ + Ai^if+^u^) + div(pV) = 0, (3.6) 

- ^{{d'' + m) ■ t^^)d'''^ - ^(d\d'' - m))d''-\ (3.7) 

p'^+^u^+i + p'^+^v^ ■ xyu"^' -fiA u^'+i + v(p'+' - /■) 
= -A(vrf'^Y(Arf'=+i - fid'^')) - Xixyd'fi^d'^' - if{d'^') - fid'))) 
- f^\u''-^ ■ + 4) - p'+V ■ yuK (3.8) 
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Define 

i^fe+lii9 II -;fc+l| 



= Wf^'Wh + P " Hi. + llvd^lh + llv^-'^^ii.. (3.9) 

Before estimates, we introduce two small positive unfixed constants 7] and e. 
From the first equation ( I3.6p . we can derive 

^llp'+^lli.<A;(t)||p'^+^||i. + r/|| V^lli^, (3.10) 

where 

Kit) = c\\ V u'm^.. + v-'c{\\ V p'mi, + mmu. 

Using the uniform estimates fl3.2l) - fl3.4l) . we obtain 

[ A'^{s)ds<C + -t, VtG[0,r,]. (3.11) 
Jo V 

fc+i 

Multiplying f l3.7p by d and integrating it over f2, we have 



— \d \ dx + \\7 d \ dx 
o-t Jn Jn 

<K(t)\\d''^'\\l^ + vi\\Vu'\\h + Wvd'Wl.), (3.12) 



where for all t G [0, T^,], 



gj(t) = Crj-'\\ V d'^+Yi^ + C\\u''-'\\l^ + Cri-'\\d^ + m\\l4d'' - m\\le 



+C||(i^ + m|Uoo||d'=-i||ioo + CT]-'\\d''-'\\l4d'' - m"2 



L6- 



The uniform estimates (I3.2p - fl3.4l) implies 

/ Bj(s)ds<C(l + -)t. (3.13) 
Jo V 



Multiplying (13. 7p by Ad , integrating it over Q and using the elliptic estimate 

l2 , we can deduce by integration by parts 

^ f \vd'^'\'dx+ I \ y^t^\Mx<CY,L^, (3.14) 
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where 



< Cr^-^\^d'^-'YH4vd'^"\\h+r^\\ V#||i., 

L2 = iim^iil6|| v^c^^^iisii v^^'^^iu^ 

< cr^-^llv'rf'^^^ll^.llv^'^'lli. + ^llv^illi., 

T II i-— 1 II II ~T^+1 II 2 

L3 = II V« llwi.ellV" IIl2) 

r II i--lll II 2~7'=+l|l II 

Li = ||m'' ^||loc|| V lU^IIVc? ||l2 

< Ce-i«^-i^.|| V rf'^'||i2 +e|| V'^'^'||i2, 



^5 = II Vc?iL3||c/'=+^ -m||Loc||d''||i6|| v^''^iL2 

< Cr'W V ^^1l^i||rf'+' -^ll^2|| V^'^'lli2 +r;|| v^l? 
= lld'^IUocll v(c^'+'-m)|U3||/|Ua|| v^'^'lU2 

< C^^-1^1l^2l| V {d'^' - Hlliill V rf'^'lli2 + ryll v^' 

L7 = IM^'IIl-IM*'^^ -m||Loo|| v^''||l2|| v^''^^I|l2 

< ^^-'ll^1l?.2|M^-^^-m||^.||v5'^'||i2 + r7||v5'lli2, 



Lg = II v(^^' + m)|U3||/+'|U6||ci'=-i|Uoc|| v5 IU2 

< c|| v(^^' + HII?.iN'"'ll//2|| v^'^'lli2, 

Lg = \\d^ + m\\HA\d^~4HA\Vd''^^\\h. 

Lio = ||/ + m||L-p''^^||L6|l Vc?''"iL3|| Vf^''^^lk2 

< C||ci'= + m||^.|| vc^'"i/fi|| V^'^'lli2, 
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Lii = WV d Wi^Wd'' - m\\Loo\\d'' ^\\Loo\\\y d ||l2 

Li2 = 0\\l4V {d'' -rn)\\L-4d''-^\\L^\\s/d''^^\\L2 

< Cv-'\\ V id' - rn)\\Ud'-4U\Vd'^'\\h+v\\vd'\\h, 
Lis = 0\\L(^\\d''^^ -m\\L^\\\y d''-^\\L:i\\\7d''^^\\L2 

< Cv-V~rn\\l4x7d'-YH4vd'^'\\h+v\\vd'\\h. 



Let's e small enough, the inequality (13 .Mp becomes 

V d \ dx + / I V " I dx 



dt 

< CC';,it)\s7t^%+Cvi\vd% + \S7u'\h), (3.15) 



where 



m\ 



c5(t) = r^\yd''+%, + \yu''-^\w^,, + \u'-%^ + r^\yd%,\d''^ 

+V'^\d%2\ V (rf'""' - m)\]j, + r^\d%,\d'''^ - m\l, 
+ 1 V {d' + m)|^i|d'^-i|H2 + + m\HAd'"%^ 
+ \d' + m|^2| V d'-%. + r'l V (t^'' - rn)\l,\d''-%2 
+r,-'\ V (^^' - + - ^Ih^I V d'-%,. 

The uniform estimates (I3.2p - (l3.4p implies 

/ Cj(s)ds <C+(C + -)t, VtG[0,r,]. (3.16) 
Jo V 

Multiplying the movement equation (13. 8p by u''^^ and integrating over Q, we can 
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deduce 



1 d 

2di 



Jn in 
Jn Jn 

+ X [ \vd''''\\Ad'+^-f(d''+^)\\u''+'\dx + X [ I v'^^ll V^''''ll^''"'|d^ 
Jn Jn 

+ X I \\/d''\\\7d''^^\\^u''+^\dx+^ I \^ dV^^ + d''\\d^^^\\d''+^\\u''+^\dx 
Jn in 

+ A / \vdV + m\\d'' -m\\t^'\\u''+^\dx+ I Ip'^+^-ZII vw'^+'ldx 
Jn Jn 
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5^M,, (3.17) 



i=l 

where we have used 



Here 



^ri\\^k+l\\2 /\\„,k-l\\2 |u,fc||2 I lU.fcll 11^^, fell ^ I ^'''11^777^+1112 

Ms <C|| V rf'^'lUHl A d'^^^ - /(^^'=+^)|U3||ll'^+i|Ue 

<^ll v«'+ii. + c|| v^'^'lli^dl v'^'+IlHI v'^'+iL6 



L2) 

y 

L2 



M4 <C 
<C 

M5 <c 
<c 

Me <C 
<C 



-jk+l. 



,2 jfell II „2 rfcii 11 „-lk+lu2 



L6, 

V^t^lUsll V^^^IIlHI Vf^''^1lL2 + 77II V^^'^'1Il2 



9' 



-fc+l||2 



V C?''IU-|| V C^^^IlHI V ^''^iLa 



-fc+1 1 

2 



A^ll ^ _fc+l||2 
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+ gllV« 



M7 <C|| Vt^1|L3||f^''+"^||L6||rf'^-m||i6||d ||i6||M^ 



Il6 



<c|| vt^lUHl v^^IIi/iN^' + '^IIhiIM'' -"^IIhiII vc^^'^^lli^ 

+ gllV« IIl2, 



M8<|b^+i-/| 



Taking Mi — Mg into fl3.17p . we obtain 
d 



V u 



k+l\2 



dx 



<CDj(t)(||/+ii. + II ^r^'u'^A\l^ + II V d'^'Wl.) + Cr^ll V 



where 



A?hA\A\w^^^ + \W^t\\LA\y u^Wl^ + v~A\^^\\lA\u^''^ 



+ II v'^^'+IlHI v'^^'+'I 



L6 



||d^+i + m||^. 



||ti'=+i-m||^i|M'=+i 



+ II c^IIlsII d^Wh^ + II V c^ll^i.6 



V d 



II V '^''llz/i 



(||d'=+^ + dA\l. Ilrf'+i^i + IM' + "^IIhi IM' - Hm) + ^'(co). 



From the uniform estimates fl3.2l) - fl3.4p . we have 

/©J(s)ds<c + (c + -)t, vtG[o,r,]. 

Jo ^ 
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Summing CTID . ([312]), f l3J[5ll and (KIM , we obtain 

Av^'^+i + (II v^-^'lli. + II v^rf'^'lli. + II v^'^^ii.^ 



< CE'^{t)^'^' + Cvi\\vu'\\i2 + \\vd Wi^), (3.20) 



where 



E;(t) = A;(t)+B;(t) + c;(t)+D;(t). 

Using the uniform estimates fl3.2l) - fl3.4l) . we obtain again 

[ E^^{s)ds <C + C{l + 7]~')t, VtG[0,T,]. (3.21) 
Jo 

Applying Gronwall's inequality to fl3.20p . we can deduce 

^'^\t) + l\\\yf^\\l. + II V^^'^'lli. + \\Vu'^'\\h)ds 
Jo 

< Cr, [ (II v^'lli^ + ll V^'lliOdrexp(C + C(l + 77-^)t). (3.22) 
Jo 

Hence choose small constants 77, T*(< T*), so that 

Cr]exp{C + C{l + r]-^)t)<^, Wt e [0,T*]. (3.23) 
We easily deduce that 

00 00 „2^* 

J2 sup ^'^\t) + J2 (llv^'^'lli. + llv^^'^'lli. + llv^'^^^lliOd. 
k=i 0^*^^* k=i Jo 

<C <oo. (3.24) 

(13.241) implies that the full sequence {p^^d^^u^) converges to a limit {p,d,u) in the 
following strong sense 

p'^ ^ p in L°°(0,T*;L2), u'' ^ u in L'^{0,T*] H^), 

-> c/ in L°°(0, T*; i/^) n L'^{0, T*; H^). 

Hence the problem (ll.ip - fll.3l) with initial boundary data fll.5p - fll.6l) has a weak solution 
{p,u,d). Furthermore using the estimates fl3.2p - fl3.4p . we obtain that a subsequence of 



27 

(p'^, , d^) converges to (p, tt, d) in an obvious weak or weak* sense. Due to the lower 
semi-continuity of various norms, from (13. 21) - (13 ■4p . {p,u,d) also satisfies the following 
regularity estimates 



sup (||p||tyi,6 + ||pt||L6 + \\u\\m + \\p\\wh^ + IbtlUo + \\d\\Hi + \\dt\\m 

Ct<T* 

+ 11 \7^u\\l2 + II \7^d\\L2 + \\\7d\\H^) 



+ I II Vpwtf + II V ut\\l2 + ||m||^2,6 + II v' dtWlz + ||d||^3dt 

< C. (3.26) 

Hence (p, w, d) is also a strong solution to the problem (ll.ip - (ll.3l) . 

4 Uniqueness and continuity in Theorem 1 

In this section, we will use energy method to prove the uniqueness and continuity 
in Theorem 1. For simplicity, we introduce some notations 

p = p — p, u = u — u, d = d — d. 

Define 

^W = llp|li^ + Plli^ + llv^lli. + llv^^lli.. 

Using the similar process in the section |3l we can obtain the following estimate (see 

(J222D) 

+ (II v^lli. + II v'^lli. + II v^lliO < cwm, (4.1) 

where F(t) G L\0,T*). 

Applying the Gronwall's inequality to (14. ip . we get for all t G [0,T*], 

m + All dWl^ + II V ^IliOdr < Cvi/(0), (4.2) 

^0 
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which imphes the uniqueness, and for all t G [0, T*], we have 

(llplli^ + IRI^. + llv^^lli.)W + f\\\ v'dWl^ + II V^lliOdr ^ (4.3) 



as {po,Uo,do) {po,uo,do) in W^'^ x H"^ x H^. 

Because d and d satisfy (11 .Sp . we obtain, similar to (13. 7p . 

^ [ I Adpdx+ [ I V^Pdx 

< c-di v^lli^ + II v^lli^ll v^iiie + II v^lliO + ||u||ioo|| a dni^. 

Applying Gronwall's inequality to the above inequality, and using the inequality (14. 3 p 
and the elliptic estimate \\d\\H2 < C\\ A d\\L2, we have, for all t G [0,T*], 

\\d\\mit)+ I II V^lli^dr^O (4.4) 
io 

as {po,uo,do) {po,uo,do) in W^'^ x H"^ x H^. 

Similarly from ( 13. 8p . using Gronwall's inequality, (I3.26P and the convergence ( 14. 3p - 
flOp . we obtain, for all t e [0,T*], 

\Mmit)+ / llVP^IIi^dt^O (4.5) 
Jo 

as {po,uo,do) {po,uo,do) in W^'^ x H"^ x H^. 

Multiplying the difference between the continuity equations by 6p^, integrating over 
(0, t) X Q and then using Gronwall's inequality, it follows from the estimate (I3.26p and 
the convergence (14.51) that for all t G [0,T*], 

||p-p|Uc(t)^0. (4.6) 

From the equations (II. 2p and (II. 3p . by a simple discussion, we can obtain for all t G 

[o,r*], 



L2 



(t), \\d\\mO,T*;H3), \\u\\l2{0,T*;H2) ^0 (4.7) 



as {po,Uo,do) {po,uo,do) in W^'^ x x H^. 

In conclusion, (I4.3p - (l4.7p complete the proof of the continuity in Theorem 1. 
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5 Proof of Theorem [2 



Suppose that there are two positive constants 9{< 1) and C such that 

max{||po||iyi>6, \\uq\\h2, ||(io - ^-11//^, hWh) < 0, (5.1) 
sup {\\v\\h^ + \\n - m\\H2 + + / || V + II^IIh/^.b 

0<t<T Jo 

+ 11 V'^tlli^ + ll^ll^adt < (7. (5.2) 

In this section, we assume the genuine constant C, maybe depending on the constant 
M(l) which is defined by fl238|) . 



By Lemma Em there exists a small 6'i(< 1) so that Vt G [0, T], V6' G (0, 6*1], 

\w^,e<C9-2, \\pt\y<Ce\, ||p|Ui,6 <C^5, |bi|U6<C^^ (5.3) 



From Lemma [2.21 we can find a small 6'2(< 1) so that G (0, ^2], 

||dt||ii(t), |M-m||^.(t), f\\d-m\\l,<\T <C^5, VtG [0,T]. (5.4) 

Jo 

By Lemma [2. 3[ a small ^3(< min{^i, 6*2}) also can be found so that \I9 G (0, ^3], 

\HH2{t), W^putWl,, [ WutWlAr, [ \\u\\lr.,edT <CeK^te[0,T]. (5.5) 

Jo Jo 

Thanks to the estimates fl5.3p - fl5.5l) . using the Theorem l2.2l we can obtain the global 
strong solution of the linear system (]2.ip -( !273|) with initial boundary value (II. 5p and 
(11. 6p provided 

max{||po||H/i.6, Ikolka, \\do - mH^s, ||fi'||i2} < 9, 

where G (0,^3]. 

Now let's talk about the iteration. First, we notice that if ^3 is taken so small that 
(76*3 < C, then the process of iteration can be continued for the same ^3. Next we will 
pay attention to the convergence of the iteration. 
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As the same process in the section 3, using the estimates f l5.3p - fl5.5l) . we can obtain, 
hke I K22^ . 

^'^\t) + All V rf'^'lli. + II v^^'^'lli. + \\vu'^'\\h)ds 

Jo 

<Cr]exp{Cit + eh + e'^ +rj-^eh + rj-^6'2)) / (II V^lli^ + II V^''lli2)dr, (5.6) 

Jo 



where 



^/c+i ^ ||-m||2^ + wt^'wh + II v^'^'lli^ + llv^^'+ii- 

Hence choose small constants r], Oq, so that V6' G (0, 6*0] and Vt G [0, T], 

Cr]exp{C{t + eh + el + r]^^eh + r]~^e^)) < ^, (5.7) 

We easily deduce that Vt G [0,T], 

00 00 „y 

J2 sup ^'^\t) + Y: / (llvrf'^'lli. + llv^rf'^'lli. + llv^^'+^lli.)d. 

0<t<T Jo 

<C<cx). (5.8) 
So we complete the proof of Theorem [2l 
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